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Abstract 

We examine the effects of the extended gauge sector of the Littlest Higgs model 
in high energy e + e~ collisions. We find that the search reach in e + e _ — > ff 
at a y/s = 500 GeV International Linear Collider covers essentially the entire 
parameter region where the Littlest Higgs model is relevant to the gauge hi- 
erarchy problem. In addition, we show that this channel provides an accurate 
determination of the fundamental model parameters, to the precision of a few 
percent, provided that the LHC measures the mass of the heavy neutral gauge 
field. Additionally, we show that the couplings of the extra gauge bosons to the 
light Higgs can be observed from the process e + e~ — > Zh for a significant region 
of the parameter space. This allows for confirmation of the structure of the can- 
celation of the Higgs mass quadratic divergence and would verify the little Higgs 
mechanism. 
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1 Introduction 



The Standard Model (SM) of particle physics is a remarkably successful theory. It provides 
a complete description of physics at currently accessible energies, and its predictions have 
been confirmed to high accuracy by all high energy experiments to date. An important piece 
of the SM remains unexplained-the mechanism of electroweak symmetry breaking. Precision 
measurements and direct searches suggest that this mechanism involves a weakly coupled 
Higgs boson with a mass in the range 114 < m# < 208 GeV at 95% CL. The Higgs mass 
parameter, however, is quadratically sensitive to UV physics. New physics at the TeV scale 
is therefore necessary to keep the Higgs light without fine-tuning. This is known as the 
hierarchy problem. Three main classes of models, supersymmetry, extra dimensions, and 
little Higgs, have been proposed to address the hierarchy problem. Which of these theories 
Nature has chosen will be determined in the coming years as the Large Hadron Collider and 
the International Linear Collider probe the TeV scale. 

The little Higgs models [TJ Ej feature the Higgs as a pseudo Nambu-Goldstone boson 
of an approximate global symmetry which is broken by a vev at a scale of a few TeV. The 
breaking is realized in such a way that the Higgs mass only receives quantum corrections at 
two loops. In contrast to supersymmetry, the one-loop contribution to the Higgs mass from 
a SM particle is canceled by a contribution from a new particle of the same spin. Little 
Higgs theories thus predict the existence of new top-like quarks, gauge bosons, and scalars 
near the TeV scale. The distinguishing features of this model are the existence of these 
new particles and their couplings to the light Higgs. Measurement of these couplings would 
verify the structure of the cancelation of the Higgs mass quadratic divergences and prove 
the existence of the little Higgs mechanism. 

The most economical little Higgs model is the so-called "Littlest Higgs" (LH) pp, which 
we introduce here and describe in more detail in Sec. 121 This scenario is based on a non- 
linear sigma model with an 577(5) global symmetry, which is broken to the subgroup 5*0(5) 
by a vev /. The vev is generated by some strongly coupled physics at a scale As ~ 47r/; 
possible UV completions of little Higgs theories are discussed in jHJH]. The 577(5) contains 
a gauged subgroup [577(2) x U(l)} 2 which is broken by the vev to the SM electroweak group 
[SU(2)l x U(1)y]- The global SU(5) breaking leaves 14 massless Goldstone bosons, four of 
which are eaten by the gauge bosons of the broken gauge groups, giving these gauge bosons 
a mass of order /. In particular, we have a heavy Z-like boson Zn and a heavy photon-like 
boson Ah which, as we will see, are phenomenologically important. The other ten Goldstone 
bosons make up a complex doublet and a complex triplet which remain massless at this stage. 
Masses for the complex triplet are generated at the TeV-scale by one-loop gauge interactions. 
The neutral component of the complex doublet plays the role of the SM Higgs. Its mass 
term comes from a Coleman- Weinberg potential and has quadratically divergent corrections 
only at two loops, giving /i 2 ~ / 2 /167r 2 . Thus the natural scale for / is around a TeV. If / 
is much higher than a few TeV, the Higgs mass must again be finely tuned and this model 
no longer addresses the hierarchy problem. 

The phenomenological implications of little Higgs models have been explored in PflEJ El 
|H1 El EH E] • Constraints arise from electroweak precision data as well as from indirect and 
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direct production at LEP-II and the Tevatron. For example, in the Littlest Higgs scenario, 
the lack of discovery of the A H , which is expected to be quite light, puts a lower bound 
on / in the few TeV range. Significant electroweak constraints come from tree-level and 
loop deviations of the p-parameter and the weak mixing angle sin 2 6 W from their SM values. 
Combining these gives a limit / > 4 TeV which is relatively parameter independent. Many 
variants of little Higgs models exist in the literature which lower this bound to / J> 1 — 2 TeV. 

In this paper we use the processes e + e~ — > f f and e + e~ — > Zh to investigate exper- 
imental limits from LEP II data on the Littlest Higgs parameters, to evaluate the extent 
of the International Linear Collider's search reach in LH parameter space, and to see how 
accurately the ILC will be able to determine the LH parameters. We will see that the ILC 
can substantially extend the discovery reach of the LHC. In addition, we will also see that 
the bounds from e + e~ — > // at LEP II exclude a large part of the LHC's search reach in 
the pp — > Zh — > Zjjn — > £ + £~bb channel. Complementary discussions of the Littlest Higgs 
model at the ILC and LHC can be found in [TUJUI]. In Sec. |2l we discuss the Littlest Higgs 
model in detail. In Sec. H3 we examine the process e + e~ — » // at LEP II and the ILC and 
determine how accurately the ILC will be able to measure the LH parameters. In Sec. 0]we 
explore the LH parameter space using the process e + e~ — > Zh at the ILC. 



2 The Littlest Higgs model and its parameters 

In this paper, we are mainly concerned with the extended neutral gauge sector present in the 
LH model. While this scenario also includes a number of parameters that arise from the top 
and scalar sectors, in which there are a number of new heavy particles, the observables of 
concern in our analysis only depend on the three parameters present in the extended heavy 
gauge sector. These are /, the vev or "pion decay constant" of the nonlinear sigma model, 
which we discussed in the Introduction, and two mixing angles. Although we focus on the 
Littlest Higgs model, we note that an enlarged gauge sector with rather generic features is 
present in all little Higgs scenarios. 

The vev / characterizes the scale of the 577(5) — ► 50(5) breaking; the effective field 
theory of the 14 Goldstone bosons has the Lagrangian 

£ s = ^Tr|P,S| 2 , (1) 

where £ is a 5 x 5 matrix parametrization of the Goldstone boson degrees of freedom QjJ ■ 
The covariant derivative contains the gauge bosons associated with the gauged subgroup 
[5*7(2) x U(l)] 2 , W x , W 2 , B u and B 2 ; 

2 

= 9 M E - i feC^-E + TWf ) + gf^B-L + ZBj)) . (2) 

At the same time, the [SU(2) x U(l)} 2 is also broken to [SU(2)l X U(1)y], and the gauge 
boson mass eigenstates after the symmetry breaking are 

W — sWi + cW 2 , W = -cW x + sW 2 , 

B = s'B x + c'77 2 , B' = -c'Bx + s'B 2 . ' ' 
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The W are the massless gauge bosons associated with the generators of SU(2)l and the B 
is the massless gauge boson associated with the generator of U{l)y- The W and B' are the 
massive gauge bosons associated with the four broken generators of [577(2) x U(l)] 2 , with 
their masses being given by 
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The mixing angles 
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, and s' = — . ^ 2 (5) 

V9JT9I V9T+9I 

relate the coupling strengths of the two copies of [SU(2) x £7(1)]. These two angles together 
with f are the three parameters of the model that are relevant to our analysis. As we 
will see, the factor of in the denominator of the expression for uib' will have important 
phenomenological consequences. 

The Higgs sector contains a scalar triplet in addition to a SM-like scalar doublet. The 
doublet and triplet both obtain vevs. The doublet vev, v, brings about electroweak symmetry 
breaking (EWSB) as in the SM, and thus v = 246 GeV. The triplet vev, v', is related to 
v by the couplings in the Coleman- Weinberg potential. Taking these to be 0(1) gives the 
relation v' ~ v 2 /2f. 

After EWSB, the mass eigenstates are obtained via mixing between the heavy (W and 
B ! ) and light (W and B) gauge bosons. They include the light (SM-like) bosons W^, Z L , 
and A L observed in experiment, and new heavy bosons W%, Zh, and Ah that could be 
observed in future experiments. At tree level, the processes e + e~ — > // and e + e~ — > Zh 
involve the exchange of only the neutral gauge bosons. Their masses are given to 0(v 2 / f 2 ) 
by 
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where m w and m z are the SM gauge boson masses, and s l 
of the weak mixing angle. Here Xh, given by [TT] 



represents the sine (cosine) 



5 , scs'c'{c 2 s' 2 + s 2 c' 2 ) 
XH ~ 2 99 5g 2 s' 2 c' 2 -g' 2 s 2 c 2 



(7) 



characterizes the mixing between B' and W' 3 in the Ah and Zh eigenstates. It is important 
to note that all but the first term in the square brackets for M\ H and M| are numerically 
insignificant. Thus M\ depends strongly on s' and not on s, and vice versa for M\ . This 
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Figure 1: (Color online) Dependence of the heavy gauge boson masses Mz H and Ma h on s 
and s', respectively, for different values of /. 



dependence is shown in Fig. [0 Note that the Ah is significantly lighter than the Zh and 
can be as light as a few hundred GeV; we will discuss the consequences of this below. 

After EWSB, the couplings of the gauge bosons Zl, Ah, and Zh to fermions similarly 
depend on s, s' and / because of the mixing between the fields. If we demand that the U(l) 
be anomaly- free, which requires y u = —2/5 and y e — 3/5 in the notation of |TT], the general 
structure of the couplings is 



g(A L ff)=g SM (Aff) 



g{Z L ff)=g SM {Zff)[l + -a l {s,s>) 



v 2 



(8) 



g(Z H ff) = ±g, 

where gsM represents the relevant coupling in the SM. A and Z are the SM photon and Z 
boson, and at and b{ are both 0(1) where i labels the species of fermion. 

The existence of the heavy gauge boson-Higgs couplings is a hallmark of the Littlest 
Higgs model. They can be probed using the process e + e~ — > Z^H through the exchange of 
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Figure 2: (Color online) The percent deviation of the vector and axial Z^ee couplings from 
the SM values for ZsM^e, taking various values for the parameters f and s. 



the Zl, Zjj, and Ah- The relevant couplings are given by 



g(Z L ^Z Lu H) = g SM (Z fl Z u H) ^1 + j^a(s, s')^ , 

g{Z Lfl Z Hv H) = — v— g^, (9) 
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where a is an 0(1) function. The formulae for the couplings can be found in Appendix B of 
the first paper in [TT] . 

Certain bounds on s and s' can be obtained by requiring that these couplings remain 
perturbative. Using the convention that a perturbative coupling g satisfies g 2 /47r < 1 gives 
s, s' > 0.1 — 0.2. Using the more conservative convention g 2 < 1 would give a smaller allowed 
range for the parameters. In the analysis that follows, we include the region where s > 0.16. 
As discussed above, expectations for the value of / arise from the requirement of naturalness. 
For / J> lOTeV, the LH model no longer addresses the hierarchy problem. 

As in [TT], we write the fermion-boson coupling as ij^igv + fiU7 5 )- It turns out that for 
the electron- coupling, \g&\ ^> while in general the shifts in the couplings due to 
mixing are roughly equal, i.e |Ag^| ~ |Agy|. Thus the relative change in gy is in general 
much greater than that for g^, as shown in Fig. El This relative change in gv is numerically 
fairly unimportant for most of the observables in our analysis, as the cross sections are 
typically functions of g v + g\. The left-right asymmetry Alh, however, has terms directly 
proportional to gv- Therefore, for the ILC, which has beam polarization capability, the A^ 
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deviation is important and introduces a surprising s' dependence in our results. We will 
discuss this in greater detail in Sec. El 

Equation El shows that for generic choices of s and s', Ma h /Mz h — s w niz / Vbmw — 1/4. 
Figure ^ illustrates this, with M^ H dipping well below 1 TeV for much of the parameter 
space. As mentioned in Sec. this light Ajj is responsible for the most stringent experimental 
constraints on the model Ej. As a result, phenomenologically viable variations of the 
Littlest Higgs models typically decouple the A H by modifying the gauge structure of the 
theory as in ^2] and I n this paper, however, we analyze the original Littlest Higgs 

model as it is the most phenomenologically well-studied. To gain some understanding of 
models in which the Ah decouples we take two approaches in our analysis. One is to choose 
a parameter value (V = a/3/5) for which the coupling of Ah to fermions vanishes. This 
decouples the Ah from all tree-level electron-positron collider physics. Another approach is 
to artificially take Ma h — ► oo while letting all other quantities in the theory take on their 
usual, parameter-dependent values. While not theoretically consistent, this approach gives 
us a more general picture of the behavior of models in which the Ah decouples. We take 
both approaches and show the results for each case throughout our analysis. 

3 Parameter Determination via e + e~^ // 

In this section we examine the process e + e~ — > ff, where all of the LH neutral gauge bosons 
participate via s-channel exchange, at past and future colliders. We first use a x- sc L uare 
analysis using the e + e~ — > // observables measured at LEP II. This analysis gives the 
region of LH parameter space excluded (to 95% confidence level) by the LEP II data. 

We then perform a similar %-square analysis at the energies and luminosity expected 
at the ILC We use the same set of observables as in the LEP II analysis as well as the 
polarization asymmetries that will be measurable due to the beam polarization capability at 
the ILC. This analysis gives the region of LH parameter space for which the ILC will be able 
to determine (to 95% confidence level) that the data cannot be explained by the Standard 
Model, and represents the ILC Littlest Higgs search reach. The two analyses just mentioned 
are described in Sec. 13.11 

Finally, in Sec. 13. 2\ we examine the ability of the ILC to determine the values of the LH 
parameters from e + e _ — > //. For a few different generic sets of LH parameters, we first 
generate sample data for the observables, and then perform a %-square analysis to map out 
the region in LH parameter space that is inconsistent (to 95% CL) with the sample data. 
The size and shape of the remaining region tells us how accurately LH parameters can be 
determined. 

3.1 The LEP II exclusion region and ILC search reach 

Here we present our numerical analysis of the experimental constraints on the Littlest Higgs 
parameter space from LEP II data as well as the search reach expected from the ILC. We 
use the Lagrangian and Feynman rules of the Littlest Higgs model as described in ^Tj. Note 
that for our analysis, we follow the notation of ^T] and take the values of the U(l) charge 
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parameters y u = —2/5 and y e = 3/5 that, as previously discussed, leave the U(l) anomaly- 
free and give the couplings shown in Eq. |S| The remaining free parameters of the model that 
are relevant to e + e~ — > ff are the sines of the two [577(2) x 2/(1)] mixing angles, s and s'; 
and the "decay constant," or vev, / as defined in Eq. [T]and Eq. |SJ 

We first study the constraints on the model from e + e _ — > // at LEP II, taking as our 
observables the normalized, binned angular distribution and total cross section for e + e~ — >bb, 
cc, and 11, with I — e, or r. We use t/s = 200 GeV and an integrated luminosity £= 627 



pb . For the detection efficiencies, we take e e = 97%, e M = 88%, e T = 49%, e b = 40%, and e c 



= 10% 14 . We perform a x-square analysis and take the SM values for all the observables 
to correspond to x 2 = 0> with a nonzero \ 2 representing deviation from the SM. This is a 
reasonable approach, since there was no detectable deviation from the SM at LEP II |T3] . 

For the ILC analysis, in addition to the above mentioned observables, we also include 
the angular binned left-right asymmetry Am for each fermion pair. We use the projected 
energy yfs = 500 GeV and luminosity C= 500 fb _1 . For the detection efficiencies, we take 
e e = 97%, e M , r = 95%, e b = 60%, and e c = 35% [T3]. 

Because of the presence of the Zh and Ah, we use a general formula for the differential 
cross section for e + e~ — > ff that is valid for any set of extra gauge bosons that can run in 
the s- channel |16j . 
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where z = cos#, Cf is the color factor, 
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and 



Bij = (viVj + l ) ( [r,r ! + , Qj = (vidj + aioAf^aj + a^Vj)^ 

Here v and a correspond to the vector and axial couplings <jy and qa discussed in Sec. 
and the sum runs over the gauge bosons in the s-channel: Al, Zl, Ah, and Zh- 

For Bhabha scattering, besides the s-channel, we also have a contribution from the t- 
channel, so that 
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Eq. ^2 i n the obvious way. 
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To calculate Alr, we need the cross sections for left- and right-handed electrons. These 
can be obtained from the above formulae by making the replacements 

1 1 

Vie — > ^ ( V ie + ^ a ie) , &ie ~> ^ ( flie + ^ ie ) ' ( 13 ) 

with A = +1(— 1) corresponding to left- (right-) handed electrons. Then the left-right asym- 
metry is given by 

A ^ z ) = r-£r&> (14) 

dz dz 

where V is the degree of e~ beam polarization at the ILC, which we take to be 80%. We 
assume the e + beam is unpolarized. 

We compute the x 2 distribution as follows, where a % represents one of the observables 
mentioned above: 

x , = Yi (^Am\\ (15) 

i ^ ' 

Here, 5a is the statistical error for each observable, given by 




(16) 



where a tat is the total cross section and ^p- is the normalized differential cross section. The 
efficiency e for each final state is given above. 

As previously noted, s,s', / are the free parameters present in the neutral gauge sector 
of the LH model. In our analysis, we choose a fixed value of s 1 and scan the parameter space 

(*,/)■ 

The exclusion region at LEP II and the search reach at the ILC correspond to the regions 
where \ 2 is greater than 5.99, representing a 95% confidence level for two free parameters. 
The combined result is shown in Fig.Elfor different values of s'. For each value of s', the LEP 
II exclusion region and the ILC search reach are on the left of the corresponding contour. 
This is because as / increases, the gauge boson masses (proportional to /) also increase (see 
Fig. H} and the deviations in the Ziff couplings (proportional to v 2 / f 2 ) decrease. For the 
ILC search reach boundary one would expect to see four contours at the upper right corner 
corresponding to the four different input values of s'. However, there is only one visible 
contour, for s' = a/3/5, because in the other three cases, the search reach covers the entire 
parameter space shown in the figure. 

As discussed above, the choice s' = a/3/5 corresponds to decoupling the A H from the 
fermion sector, as verified by the results shown in Eq. |S| In this case, the ILC search reach 
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95% CL contours with light A H at Vs = 500 GeV 
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Figure 3: (Color online) LEP II exclusion region from e + e~ — * // and ILC search reach in 
the parameter space (s, /) for different input values of s', and including the Ah contribution. 
For s' = a/3/5 there are two lines of the same symbol/color, one on the boundary of the 
LEP II exclusion region, and one on the boundary of the ILC search reach region. For the 
other values of s' the curve shown is the boundary of the LEP II exclusion region, while the 
ILC search reach covers the entire parameter region shown. 



can be as low as / ~ 2 TeV for large values of s. For other values of s', the search reach 
is greater than / ~ 10 TeV for all values of s. We thus see how strongly the presence of 
the relatively-light Ah can affect the phenomenology. For LEP II, the story is similar; the 
exclusion region for s' = a/3/5 is much smaller than for other values of s'. This is because 
the observed deviation at s' = a/3/5 is solely due to the modification in the Z L ff coupling 
and the presence of the Zh, which is generally several times heavier than the Ah- For other 
values of s' the constraints on / can be as high as ~ 6 TeV. Overall, the LEP II exclusion 
regions have constraints on the parameter / that are roughly consistent with those from 
precision measurements [TJ El 13 El E] • 

As discussed in Sec. 121 we also examine the general behavior of models in which the Ah is 
decoupled by taking Ma h —>■ oo while letting all other quantities take on their usual values. 
The results in this case are presented in Fig. 0] It is not surprising that the s' = a/3/5 
contours in Fig. HI are exactly the same as in Fig. El since the Ah is decoupled for this 
choice of s'. For other values of s', the contours are very different in the two cases. The s 
dependence of the contours in Fig. HJis easy to understand. The Z H ff couplings go as gc/s, 
thus the ILC contours show that the search reach is higher for lower values of s. Similarly, 
for LEP II, the exclusion region extends farther out in / for lower values of s. There is, 
however, a "turnover" for the LEP II exclusion region around s ~ 0.3 where the contours 
start moving towards lower /. This takes place because the mass Mz H begins to increase 
(see Fig. [TJ) and the overall contribution from Z H to the observables starts to decrease as s 
gets smaller. 
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95% CL contours with M A -> °o at Vs = 500 GeV 
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Figure 4: (Color online) LEP II exclusion region and ILC search reach as in Fig. El but with 
M Ah -> oo . 



With Ma h — > oo, the s'-dependence of the \ 2 is mostly due to the deviation in the Z^f f 
couplings, since neither the Z H f f couplings (see Eq. EJ) nor M Z[1 (see Eq. EJ) are strongly 
dependent on s'. This explains why there is less variation among the different contours in 
Fig. 0] than in Fig. El For values of s close to 1, however, the ILC contours for different 
values of s' begin to deviate from one another markedly. This s' dependence is due to the 
s'-sensitive deviation of Alr, as discussed in Sec. El This is confirmed by Fig.0 which shows 
the relative contribution of the different observables to the \ 2 & t the ILC with Ma h — > oo 
and s = 0.95. Note that Am for various final states dominates the \ 2 where it is large. 

The fact that the search reach is lowest for s' = a/3/5 then indicates that the deviations 
in the Z^ff couplings are smallest for this parameter value. This coincidence arises because 
both An and Z L are linear combinations of gauge eigenstates. An to lowest order is just 
B', whose couplings to fermions vanish at s' = a/3/5. As the s'-dependent terms in the 
deviation of the Z^ff couplings arise from the B' admixture, they also vanish at this value. 
This is also confirmed by Fig. which shows that the relative contribution of A^r and the 
total x 2 decrease around s' = a/3/5. 

The search reach at a a/s = 1 TeV ILC reaches to somewhat higher values of the pa- 
rameter s, but has essentially the same reach for the parameter / as the a/s = 500 GeV 
machine. This is reasonable; as s approaches unity, the contribution from the deviations in 
the Zl couplings dominates the search reach, and these contributions are not as important 
as the center of mass energy increases. The result is that the search reach is very simlar for 
both y/s = 500 GeV and 1 TeV. We will see later, however, that the \fs = 1 TeV data can 
significantly improve the parameter determination. 

Figure El and Fig. |U show that the yfs = 500 GeV ILC search reach in general covers 
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Contribution to % 2 (f=6.0 TeV, s = 0.95, M Ah ^ ~) 
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s' 

Figure 5: (Color online) Fractional contribution to the total x 2 f° r each e + e~ — > ff ob- 
servable at a 500 GeV ILC for fixed (s, /). The labels on the legend go from top down in 
one-to-one correspondence with the shaded sections. For example, Alr(c + e~) contributes 
~ 78% to the total x 2 at s' = 0.4. The line labeled "95%" is the total x 2 /5.99. This means 
that the region s' ~ [0.55, 0.9] where this line dips below 1 is outside the y/s = 500 GeV ILC 
search reach. 



most of the interesting parameter space where the Littlest Higgs models are relevant to the 
gauge hierarchy. Thus the e + e~ —>■ ff process alone is an effective tool for investigating the 
Littlest Higgs model at a ILC. 

It is important to compare the ILC search reach to that of the LHC. An ATLAS based 
analysis of the LHC search reach for the heavy gauge boson Zh of the Littlest Higgs model 
was computed in Ref. [T7]. The resulting 5a contour for discovery of the Zh at the LHC is 
reproduced in Fig. |H1 (where we have converted their results to our choices of parameters / 
and s for the axes). Fig. also displays our results for the ILC (taking Ma h — > oo), where 
we have now employed a statistical significance of 5cr rather than 95% to facilitate an equal 
comparison. We see that the ILC substantially extends the LHC search reach for s ^ 0.8. 

3.2 Parameter Determination: sample fits 

We have now determined the available parameter space accessible to the ILC and not already 
excluded by LEP II. It remains to ask, given the existence of an LH model with parameters 
in this accessible range, how accurately would the ILC be able to measure them? It is well- 
known [THJ HH] that the ability to precisely measure the couplings of heavy gauge bosons is 
one of the fortes of the ILC. 

We first discuss the capability of the LHC to determine the LH model parameters. Nu- 
merous studies jIHl EOj have addressed the ability of the LHC to determine the couplings of 
new gauge bosons. The results of these studies show that while some model differentiation is 
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Figure 6: (Color online) Comparison of ILC and LHC search reach. The LHC data was 
taken from Fig. 8 of The search reach lies to the left of and underneath the contours. 



possible for Z' bosons with masses % 2 TeV, absolute determination of the Z' couplings is 
not possible. There are three main reasons for this: (i) the number of observables is limited 
in the hadron collider environment. The observables are the number of events (i.e., cross 
section times branching fraction), the forward-backward asymmetry, and the rapidity asym- 
metry) for leptonic final states only. Other final states are not detectable above background, 
(tt final states are a possible exception, but such events will be heavily smeared and thus 
not useful for a coupling analysis), (ii) The observables are convoluted with all contributing 
parton densities, (iii) The statistics are insufficient for Mz< ^ 2 TeV. Here, in the case of 
the LH, our results show that LEP II essentially excludes the region Mz H ^ 2 TeV, and 
thus we do not expect the LHC to contribute to the parameter extraction in any significant 
way. We note, however, that a very precise mass measurement for Zh will be obtained at 
the LHC. 

To determine the accuracy of the parameter measurements, we perform some sample fits, 
using a x-square analysis similar to the one described in the preceding section. We use the 
same ILC observables as before. In some cases we also include data from a y/s = 1 TeV 
run, for which we also take an integrated luminosity of £= 500 fb~ . We note that we can 
exchange Mz H for /, so we now take Mz H , s, and s' as our free parameters. We choose 
a generic data point (s, s' , M Zh ) and use it to calculate the observables, which we then 
fluctuate according to statistical error. We assume that the Large Hadron Collider would 
have determined M Z[I relatively well, to the order of a few percent for M Zh ^5 — 6 TeV; we 
thus fix Mz H and perform a 2-variable fit to s and s' . Scanning the s-s' parameter space, 
we calculate the \ 2 a ^ every point. We find the minimum \ 2 point; the 95% CL region 
surrounding it is the region for which the x 2 is l ess than this minimum \ 2 plus 5.99. 

Figure [7| shows the results of this fit for two sample data points in the contrived scenario 
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Sample fits for 95% CL (M Zh = 3.0 TeV, °°, Vs = 500 GeV) 



Data point • 

s'=V3/5 

s'= 0.5 
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s' 

Figure 7: (Color online) 95% CL sample fits to the data points (s = 0.5, s' = 0.5) and (s = 
0.5, s' = y/3/5), using e + e — > // observables at a 500 GeV ILC, taking Mz H = 3.0 TeV 
and Ma h — > oo. 



with Ma h — > oo. For both of these points, the determination of s is very accurate. This 
is due to the strong dependence of the Z H f f couplings on s, as discussed in the previous 
section. The s' determination is worse than that for s because of our choice s = 0.5. At this 
value of s, the contributions from the Zl coupling deviations (which carry the s' dependence) 
are smaller than the Zh contributions. The reason the s' determination is better for s' = 0.5 
than it is for s' = a/3/5 is that the s'-dependent deviations in g Vz f j vanish for the latter 
value. 

Figure |Hl shows the results from a similar fit and illustrates how it can be improved with 
data from a higher-energy run with yfs = 1 TeV at the ILC. Here, the s determination is 
not much more accurate than the s' determination, as the s'-independent Zh contributions 
no longer dominate the fit for s = 0.9. 

In Fig. El we show results from a fit with the full Ah contributions. Not surprisingly, 
the parameter determination is much more precise, as the Ah contributions, when present, 
dominate the x 2 . Since the Ah couplings depend only on s', it is also no surprise that here 
the s' determination is in general much better than that for s. 

If, for some reason, the LHC doesn't provide a good measurement of M Zh , we would 
need to include that quantity, or equivalently /, in our fits to the data. In Fig. EH we show 
the results where we have set s' = a/3/5 and we fit to s and /. Note that for one of the 
data points, the allowed region doesn't close. This highlights the importance of using both 
the LHC and the ILC to reliably determine the model parameters. 
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Sample fits for 95% CL (M Zh = 3.3 TeV, M A °°) 
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Figure 8: (Color online) Like Fig. 0except Mz H = 3.3 TeV and the data points are (s = 0.9, 
s' = 0.5) and (s = 0.9, s' = a/3/5). Also shown for each point is an improved fit from adding 
data from a y/s = 1 TeV, £= 500 fb _1 run at the ILC. 



Sample fits for 95% CL (M Zh = 3.3 TeV, Vs = 500 GeV, light A H ) 
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Figure 9: (Color online) Like Fig. except M Zh = 3.3 TeV and the data points are (s = 0.65, 
s' = 0.65) and (s = 0.65, s' = a/3/5), and the full M& H contributions are included. 
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Sample fits for s'=V3/5 
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Figure 10: (Color online) Sample fits to the data points (s = 0.5, / = 4 TeV) and (s = 0.9, 
/ = 4 TeV), taking s' = ^/3/5. At the decoupling value of s', the A H does not contribute. 

4 Parameter determination using e + e~ — > Z^/i 

In order to confirm that the LH model is the correct description of TeV-scale physics, it is 
important to test the hallmark of the LH mechanism, namely that the Higgs mass quadratic 
divergences are canceled by new particles with the same spin as their SM counterparts. 
The proof lies in the measurement of the new particle couplings to the Higgs. Here we are 
concerned with the coupling of the heavy Z to the Higgs boson. This coupling can be tested 
via the process e + e~ — > Z^h. In the LH model, deviations of observables related to this 
process from their SM expectations come from three sources: the diagram with the Zh in 
the s-channel, the diagram with the Ah in the s-channel, and the deviation of the Z^Z^h 
coupling from its SM value. 

In this section we repeat the analysis of Section using the process e + e~ — ► Z^h and 
taking the total cross section as our observable with mh = 120 GeV. We assume that at a 
y/s = 500 GeV ILC this cross section will be measured to an accuracy of 1.5% [T5j . 

The cross section, including the effects of additional gauge bosons, can be written as 



°z L h = 




where P?? was defined in Eq. El The sum runs over the participating gauge bosons in the 
s-channel: Zl, Ah, and Zh- Here, |k| is the magnitude of the 3-momentum of the outgoing 
Z L , given by 

|k| = ^j= s s]{m 2 H ~ M\ L Y + a(a - 2(Mf L + ml)). (18) 




\ 

\ Input point • 
\ s=0.9,Vs = 0.5 TeV 

\ s=0.5,Vs = 0.5 TeV 

s=<\9,Vs = 0.5,1 TeV (combined) 
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Figure 11: (Color online) LEP II exclusion region from e + e~ — ► // and ILC search reach in 
the parameter space (s, /) from the process e + e~ — > Z L h, for various values of s' and includ- 
ing the full Ah contributions. For each value of s' there are three curves; one corresponds 
to the LEP II exclusion region, and the other two represent the ILC search reach region for 
y/s = 500 and 1000 GeV, respectively, taking an integrated luminosity of 500 fb _1 at each 
center-of-mass energy. 



The couplings V{ and a\ are the same as before-the axial and vector couplings of electrons 
to the zth gauge boson. 

We carry out the x-square analysis as before. Figure ^2 shows the ILC search reach 
in the LH parameter space, where each plot corresponds to a different choice of s'. By 
comparing to Fig. H] we note that e + e~ — > Zjjn gives a much poorer search reach than 
e + e~ — > //. In particular, when s' is near the decoupling value \/3/5 the LH model is 
generally indistinguishable from the SM. Well away from s' = a/3/5, as shown for s' = s/2 
and s' = 0.5, the search reach covers almost all of parameter space, except for regions of low 
/ where interference between the Ah and Zh conspire to bring the cross section near its SM 
value. These regions, however, are ruled out by LEP. 

In the case s' = s, however, there are regions that exhibit similar interference effects and 
are not ruled out by LEP data. For example, consider the two data points / = 4.0 TeV, 
s = 0.61 with (a) s' = and (b) s' = 0.61. With y/s = 500 GeV, (b) is within the 

search reach while (a) is just outside the search reach. Figure shows that at this value 
of y/s, the deviation of the cross section from the SM is much greater for s' = 0.61 than for 
s' = a/3/5, since the A H decouples in the latter case. With y/s = 1 TeV, this behavior is 
reversed; point (a) is outside the search reach while (b) is within. At this value of y/s the 
interference between Ah and Zh brings the cross section closer to the SM value when the 
Ah contributes. 

Figure EH shows the search reach obtainable with 500 fb _1 at a y/s = 500 GeV ILC, taking 
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Figure 12: (Color online) The cross section for e + e~ — > Zih as a function of y/s for the 
SM and two different points in LH parameter space. The insets show the behavior near the 
expected ILC y/s values of 500 GeV and 1 TeV. The resonance at about 700 GeV corresponds 
to the Ah- 



Ma h — > oo . Comparing to Fig. El we see that the search reach here is much smaller than for 
e + e~ — > //. Figure ITU displays the corresponding reach at y/s = 1 TeV with 500 fb _1 . In 
both cases, and for all choices of s', the search reach decreases markedly around s = l/y/2. 
This is because the Z L Z H H coupling vanishes at this value of s, as can be seen in Eq. 
It is also interesting to note that the spread in the search reach as s' is varied is larger for 
y/s = 500 GeV than it is for 1 TeV. This can be understood if one notes that y/s = 1 TeV 
is closer to the Z H pole (as M Zh — a few TeV throughout the parameter space) than 
is 500 GeV. Thus the deviation of o Zh from its SM value at y/s = 1 TeV is dominated 
by the presence of the Z H , whose mass and couplings are essentially s'-independent. At 
y/s = 500 GeV, the deviation of ozk has a more significant contribution from the deviation 
of the e + e~Z L coupling, which is strongly dependent on s' (see Fig- EJ) - For both values of 
y/s, the sensitivity in the range of parameter space where s ^ 0.5 does not reach beyond the 
general precision electroweak bound of / J> 4 TeV. 

One could hope to improve the sensitivity by adding the measurement of the Higgs 
branching ratios as additional observables. It turns out, however, that the LH deviations 
of the branching ratios from their SM values are at most 1-2%, which is smaller than or 
equivalent to the experimental sensitivity at the ILC. 

Lastly, we again compare the reach obtainable at the ILC from this process to that of 
the LHC in pp ZhZl + h. We display the 5a results from the ATLAS based analysis 
[T7j of this process in the LH using the final state £ + £~bb in Fig. E3 We also show our 
results, again adjusted for 5cr rather than 95% statistical significance. This figure shows 
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95% CL contours with M A -> °° at Vs = 500 GeV 
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Figure 13: (Color online) The ILC search reach from the process e + e~ — » Z^h for various 
values of s', taking ^/s = 500 and Ma h — > oo. The LEP II exclusion region from e + e~ — > // 
is shown for s' = s/2 (from Fig. HJ) for comparison. 
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Figure 14: (Color online) Same as Fig. EH but for y/s = 1 TeV. 
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Figure 15: (Color online) Comparison of ILC and LHC search reach in the Zh — ► Zh channel. 
The LHC curve was read from Fig. 22 of |17j . 



that the ILC overwhelms the capability of the LHC in this channel. In fact, our analysis of 
e + e~ —>■ ff shows that for s ^ 0.8 the LEP II results already exclude the possibility of the 
LHC observing the Z L + h decay of the Z H . 

5 Summary 

Little Higgs models provide an interesting mechanism for addressing the hierarchy problem. 
They contain a single light Higgs boson which is a pseudo-Goldstone boson with a small 
mass generated at the two-loop level. The quadratically divergent loop contributions to the 
mass of this Higgs are canceled by contributions from new particles appearing at the TeV 
scale. These cancelations take place between contributions from particles which have the 
same spin. Measurement of the couplings of these new particles to the light Higgs would 
verify the structure of these cancelations and establish the Little Higgs mechanism. 

Here, we have investigated the extended gauge boson sector within these theories. Numer- 
ous Little Higgs models, based on various global symmetries, have been proposed. However, 
the existence of an enlarged gauge sector, with rather generic features, is endemic to all 
these scenarios. We choose to work in the framework of the simplest model of this type, 
known as the Littlest Higgs, based on an SU(5)/SO(5) nonlinear sigma model. This scenario 
contains the new heavy gauge bosons W H , Zh, and Ah in addition to the SM gauge fields. 
The masses of these additional gauge bosons are expected to be of order the global symme- 
try breaking scale of / ~ TeV. (It is expected that / < 10 TeV in order for this scenario 
to be relevant to the hierarchy.) However, due to the group theory structure, the Ah can 
be significantly lighter resulting in stringent constraints from precision electroweak data. 
Phenomenologically viable Littlest Higgs models must thus decouple the Ah and we have 
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examined two such approaches in our analysis. One, where we choose the model parameters 
such that the fermion couplings of the Ah vanish, and another where we artificially take 
M Ah -> oo. 

We study the effects of the new neutral gauge bosons in e + e~ annihilation. These particles 
can participate in e + e~ — > ff and e + e~ —>■ Zh via s-channel exchange, and their effects can 
be felt indirectly for center of mass energies well below their masses. We find that fermion 
pair production is more sensitive to Little Higgs effects than Zh associated production. We 
perform a thorough investigation of the model parameter space and find that observables at 
LEP II exclude the region / < 1 — 3 TeV, which is consistent with the constraints obtained 
from precision electroweak data. The search reach of the proposed International Linear 
Collider, operating at ^fs = 500 GeV, covers essentially the entire parameter region where 
this model is relevant to the hierarchy, i.e., f ^ 6 — 10 TeV. In the case of a 1 TeV ILC, the 
search region probes slightly larger values of the mixing parameter s, but similar values of 

We have also demonstrated that once a signal is observed in these channels, accurate 
measurements of the couplings of the heavy gauge fields can be obtained from fermion pair 
production at the ILC. These couplings are related to the mixing angles in the extended 
gauge sector and we show that experiments at the ILC can determine the fundamental 
parameters of the theory. For illustration, we performed a fit to generated data for sample 
points in the Littlest Higgs parameter space, and found that the fundamental parameters can 
be determined to the precision of a few percent, provided that the LHC measures the mass 
of the heavy neutral gauge field. If information on the new boson masses is not available 
from the LHC, then the parameter determination at the ILC deteriorates. Additionally, the 
couplings of the extra gauge bosons to the light Higgs can separately be determined from 
e + e~ — > Zh for a significant region of the parameter space. This enables ILC experiments to 
test the consistency of the theory and verify the structure of the Higgs quadratic divergence 
cancelations. 

In summary, we find that the ILC has the capability to discover the effects of the Littlest 
Higgs model over the entire theoretically interesting range of parameters, and to additionally 
determine the couplings of the heavy gauge bosons to the precision of a few percent. 
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